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Abstract 

An  engineering  model  that  can  be  used  to  obtain 
predictions  of  axial  distributions  of  temperature  and 
species  concentrations  in  complex  flows  has  been  for¬ 
mulated  and  applied  to  waste  incineration  in  a  dump 
combustor.  The  model  incorporates  mean  convec¬ 
tion  and  molecular  diffusion  in  a  quasi-one  dimen¬ 
sional  sense  and  uses  a  stochastic  model  to  approxi¬ 
mate  the  effects  of  turbulent  convection.  The  inputs 
to  the  model  are  extracted  from  experimental  data 
and  results  of  Large  Eddy  Simulations  (LES).  Com¬ 
parisons  of  waste  consumption  rates  and  pollutant 
formation  in  the  dump  combustor  with  experimental 
data  show  that  the  model  can  capture  the  correct 
trends  and  achieve  fairly  good  quantitative  agree¬ 
ment.  Further,  the  model  does  not  use  any  ad-hoc 
constants  that  need  to  be  calibrated  or  tuned. 

1  Introduction 

There  is  a  strong  need  in  the  industry  for  meth¬ 
ods  that  accurately  analyze  combustor  flows  in  order 
to  obtain  quick  solutions  to  design  problems.  Sim¬ 
ple  models  or  methods  capable  of  predicting  mean 
values  without  undue  computational  expense  are 
needed.  An  additional  feature  desired  of  these  mod¬ 
els  is  that  they  should  be  capable  of  handling  var¬ 
ied  flow  conditions  without  requiring  recalibration. 
These  models,  when  developed,  can  then  be  used  to 
parameterize  the  flow  in  complex  geometries  with  re¬ 
spect  to  any  of  the  variables  of  the  problem.  In  this 
paper,  we  address  this  need  by  developing  a  simple 
model  to  predict  complex  combustor  flows. 

Many  of  the  earlier  attempts  to  model  com- 
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plex  geometry  combustor  flows  using  simplified  mod¬ 
els  have  been  largely  directed  towards  the  control  of 
instabilities  that  develop  in  these  combustors[l,  2,  3]. 
Most  of  these  methods  use  some  procedure  that  sim¬ 
plify  the  equations  (typically,  linearization)  to  ob¬ 
tain  equations  for  the  perturbations  in  the  flows. 
More  often  than  not,  this  assumes  knowledge  of  the 
mean  flow  field  information  to  solve  the  perturbation 
equations.  Yang  and  Culick [4]  used  an  integral  tech¬ 
nique  to  solve  for  the  mean  flow  variables  and  then 
used  this  information  to  solve  the  linearized  pertur¬ 
bation  equations.  However,  this  analysis  was  used 
to  mainly  study  the  modes  of  acoustics  and  insta¬ 
bilities  in  a  liquid  fueled  ramjet.  In  another  work, 
Logan  et  al.  [5]  analyzed  the  acoustic  modes  of  a  low 
speed  dump  combustor.  Again,  the  mean  flow  field 
information  was  assumed  known  rather  than  calcu¬ 
lated.  In  an  earlier  work,  Jou  and  Menon  [2]  used 
the  information  obtained  from  numerical  simulations 
to  construct  a  model  for  the  pressure  oscillations  in 
a  ramjet.  ,  Sterling  et  al.  [6]  studied  the  longitu¬ 
dinal  mode  instabilities  in  a  dump  combustor  using 
experimental  measurements.  Thus,  although  signif¬ 
icant  amount  of  work  has  been  directed  towards  the 
phenomena  of  instabilities  and  pressure  oscillations 
in  combustor  flows,  relatively  little  effort  has  gone 
into  developing  simple  models  to  predict  the  mean 
reacting  flow  fields  in  these  combustors. 

In  this  paper,  results  from  an  ongoing  re¬ 
search  initiative  to  develop  an  engineering  model  to 
study  complex  combustor  flows  are  presented.  The 
model  is  being  developed  to  study  the  controlled 
incineration  of  toxic  wastes  in  a  dump  combustor. 
This  approach  parallels  an  ongoing  experimental 
study  of  an  identical  dump  combustor  [19].  The 
basic  model  formulation  and  details  are  presented 
in  the  next  section.  The  model  uses  the  concept 
of  stochastic  stirring  events  to  model  the  effects  of 
turbulent  stirring,  first  postulated  by  Kerstein  [12], 
in  conjunction  with  a  quasi-one  dimensional  solver. 
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The  model  uses  as  input,  parameters  that  describe 
the  turbulence  in  the  flow  field  such  as,  the  inte¬ 
gral  length  scale  variation  and  the  rms  of  the  ve¬ 
locity  fluctuations.  These  parameters  are  obtained 
from  experimental  measurements  of  the  turbulent 
flow  field.  The  parameters  governing  the  stirring 
events  are  derived  based  on  three  dimensional  scal¬ 
ing  laws  for  homogeneous  turbulence. 

2  Model  Formulation 

In  the  approach  adopted  here,  the  turbulent  flow 
field  in  a  dump  combustor  is  conceptually  broken 
down  into  the  various  physical  processes  occurring 
in  the  combustor  shear  layer.  The  processes  con¬ 
sidered  are  the  following  :  1)  Mean  convection  at 
the  convection  velocity  of  the  structures  in  the  shear 
layer,  2)  The  turbulent  convection/stirring  due  to 
the  presence  of  vortices  in  the  shear  layer,  3)  Molec¬ 
ular  diffusion  due  to  the  presence  of  scalar  gradients, 
and  4)  chemical  reaction  and  the  accompanying  heat 
release.  Each  one  of  these  physical  processes  is  han¬ 
dled  distinctly  and  these  various  models  and  their 
relevant  parameters  are  discussed  below. 

The  mean  convection  of  the  scalars  in  the 
shear  layer  is  handled  using  the  quasi-one  dimen¬ 
sional  equations.  The  reduced  convection  velocities 
in  the  shear  layer  are  modeled  as  the  effects  of  area 
changes.  It  should  be  noted  here  that  we  are  in¬ 
terested  only  in  modeling  the  turbulent  mean  flow 
and  not  the  unsteady  flow  variations.  Thus,  by  us¬ 
ing  either  full  numerical  simulation  or  experimental 
measurements  (as  is  done  here)  the  target  velocity 
profiles  are  known  apriori.  The  area  distribution 
is  adjusted  until  the  mean  flow  velocity  profile  is 
obtained.  Large  variations  are  not  expected  in  the 
mean  velocity  profiles  for  small  changes  in  the  oper¬ 
ating  parameters  of  the  system.  Hence,  this  calibra¬ 
tion  is  expected  to  be  required  only  once  for  a  given 
geometry. 

The  second  aspect  concerns  modeling  of 
fuel-air  mixing.  Fuel  injection  is  modeled  as  simple 
mass  addition  into  the  variable  area  flow.'  This  gives 
rise  to  source  terms  in  the  governing  equations  for 
mass,  momentum,  energy,  and  mass  fractions  of  the 
species.  Mixing  is  modeled  in  two  parts.  Molecular 
mixing  is  accounted  for  by  directly  including  the  dif¬ 
fusive  terms  in  the  governing  equations  for  the  quasi- 
one  dimensional  flow.  Fickian  diffusion  is  assumed 
for  the  species  and,  in  this  work,  all  the  species  are 
assumed  to  have  equal  diffusivities  and  unity  Lewis 
number  although  the  formulation  is  generic  . 

Turbulent  stirring  is  accounted  for  by  the 


stochastic  stirring  events  that  mimic  the  action  of 
an  eddy  on  a  scalar  field.  This  is  done  by  random 
rearrangement  events  of  the  scalar  field  along  the 
domain.  The  laws  governing  these  rearrangements 
are  derived  such  that  the  apparent  diffusion  induced 
by  these  processes  mimic  that  experienced  by  a  par¬ 
ticle  in  a  real  turbulent  flow.  In  order  to  do  this,  the 
high  Reynolds  number  scalings  based  on  the  Kol¬ 
mogorov  cascade  are  used.  The  stirring  events  can 
then  be  described  in  terms  of  two  parameters.  The 
first  one  is  a  stirring  frequency  parameter  (A)  which 
determines  the  rate  at  which  these  stirring  events 
occur.  The  second  parameter  is  a  pdf  describing  the 
size  distribution  of  the  rearrangement  events.  The 
actual  expressions  for  these  parameters  depend  on 
the  kind  of  mapping  event  chosen.  In  this  study  the 
block  inversion  mapping  is  chosen.  Mathematically, 
the  block  inversion  mapping  of  an  arbitrary  scalar 
field  9(x1t)  can  be  expressed  as  follows, 

fit*  4\  _  /  ~  if  /©  —  /  <  x  <  /0  -|-  / 

'  ’  ~  \  9(x,t)  otherwise 

(!) 

Detailed  description  of  the  mapping  and  its 
diffusive  and  dispersive  properties  are  discussed  by 
Iverstein  [12]  and  are  omitted  here  for  brevity.  Using 
the  derivations  given  there,  the  following  expressions 
are  derived, 

24l/fle,Ll  -(!/»?)— 5/3 
5  7 f  (i/77)4/3  -  1  V  1 

c  r-  8/3 

^=3^/3-L-/3  <8> 

where,  v  is  the  kinematic  viscosity,  L  is  the 
integral  length  scale,  Rej  is  the  Reynolds  number 
based  on  the  rms  velocity  fluctuations  and  the  in¬ 
tegral  length  scale,  and  77  is  the  kolmogorov  length 
scale.  Note  that,  in  general,  since  L  and  urms  are 
also  a  functions  of  space,  both  A  and  /(/)  are  in  gen¬ 
eral  functions  of  space. 

The  quasi-one  dimensional  equations  are 
marched  in  time.  When  the  epoch  for  a  stirring 
event  occurs,  the  location  of  the  event  is  chosen.  In 
order  to  do  this  we  define  the  global  event  rate  E(t) 
at  time  t  as 

m  =  j  A  (x,t)dx  (4) 

Now,  the  event  location  is  selected  randomly 
with  probability  \(x,t)dx/ E(t)  that  its  center  falls 
in  the  interval  (x,x  +  dx).  Once  the  location  is  se¬ 
lected,  the  event  size  is  chosen  by  randomly  sampling 
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the  pdf  /(/).  Then,  the  scalar  field  is  subjected  to 
the  mapping  given  in  1.  This  is  done  in  such  a  way 
that  equal  masses  are  interchanged  between  cells  un¬ 
dergoing  the  inversion  process.  After  the  inversion 
at  the  current  time  has  been  completed,  the  time 
interval  to  the  next  event  is  determined  by  sam¬ 
pling  an  exponential  distribution  with  mean  1/E(t). 
The  solver  then  goes  back  to  the  quasi-one  dimen¬ 
sional  system  until  that  time  for  the  next  inversion 
is  reached  at  which  point  the  above  process  is  again 
repeated.  The  aspects  of  this  implementation  spe¬ 
cific  to  its  application  to  the  shear  layer  in  a  dump 
combustor  are  discussed  later. 

The  final  aspect  is  the  modeling  of  the  chem¬ 
istry  in  the  combustor.  In  the  cases  considered  here, 
a  mixture  of  ethylene  and  benzene  is  used  as  fuel  and 
waste.  Ethylene  is  the  actual  fuel,  while  the  benzene 
is  used  as  a  waste  surrogate.  This  study  is  a  part  of  a 
larger  effort  to  evaluate  the  use  of  a  dump  combustor 
type  geometry  for  the  purpose  of  waste  incineration. 
In  order  to  evaluate  fuel  and  waste  consumption 
rates  in  the  combustor,  a  realistic  reaction  mecha¬ 
nism  has  to  be  used.  The  complete  mechanism  for 
the  combustion  of  ethylene  consists  of  277  reactions 
involving  48  species[7].  It  is  computationally  forbid¬ 
ding  to  use  this  full  mechanism.  Hence,  a  methodol¬ 
ogy  suggested  by  Singh  and  Jachimowski[8]  is  used. 
This  reduces  the  mechanism  to  a  more  manageable 
11  reactions  among  10  species.  Pollutant  formation 
in  combustors  of  the  type  considered  here  is  a  major 
problem  that  needs  to  be  addressed.  Therefore  an 
additional  reaction  based  on  the  suggestion  by  Chen 
and  Kollmann[9]  for  the  formation  of  NOx  is  used  for 
this  purpose.  The  final  reaction  system  contains  12 
species  and  12  reactions  and  is  given  in  Table  l.The 
coefficients  from  CHEMKIN  [11]  chemical  kinetics 
package  is  used  to  calculate  the  free  energies  for  the 
computation  of  the  reaction  rate  data. 

In  the  following  section,  the  equations  used 
in  the  model  are  presented.  The  various  sub-models 
described  above  are  discussed  in  greater  detail  and 
methods  of  their  calibration  are  discussed. 


3  Governing  Equations 

The  governing  equations  are  the  basic  conservation 
laws  for  mass,  momentum,  energy  and  species  mass 
fractions  for  a  quasi-one  dimensional  flow  with  area 
variation.  These  equations  for  a  multi-component 
system  are  written  as  follows: 


dQA 

dt 


+  ^A=  Ga  +  Gm 

OX 


(5) 


Reaction 

A 

n 

Ea 

C2H4  -j-  3O2  200  4"  2H2 

2.2ell 

0.0 

1.5558e8 

C6  H6  +  3 O2  6 CO  +  3H2 

1.8e8 

0.0 

1.4859e8 

CO  +  0  5±  C02  +  M 

5.3el 

0.0 

-1.9004e7 

CO  +  OH  C02  4*  H 

4.4 

1.5 

-3.0975e6 

H2  4-  02  ^  OH  4-  OH 

1.7e7 

0.0 

2.0092e8 

H  +  Oi^tOH  +  O 

2.6e8 

0.0 

7.0322e7 

OH  4-  Hi  HiO  4-  H 

2.2e7 

0.0 

2.1557e7 

0  +  Hi?±0H  +  H 

1.8e4 

0.0 

3.7254e7 

OH  4-  OH  ^  HiO  4-  0 

6.3e7 

0.0 

4.5626e6 

H  4-  H  ^  Hi  4-  M 

6.4e5 

-1.0 

0.0 

H  +  OH^HiO  +  M 

2.2el0 

-2.0 

0.0 

Ni  +  Oi^  2 NO 

1.82e8 

0.0 

1.6061e8 

Table  1:  Chemical  reaction  mechanism  used  in  the 
present  work  and  in  the  LES  simulations 


pu 

pu2+p 

puH 

pYkU 


(6) 

(7) 

(8) 


(9) 


In  the  above,  p  is  the  density,  u  is  the  ve¬ 
locity,  E  is  the  total  energy  per  unit  mass  given  by, 


(10) 


where  e  is  the  internal  energy  per  unit  mass  given 


by, 

e  =  h-P- 
P 

(11) 

and 

N. 

h  =  YJykhk 

(12) 

k= 1 

hk  =  Ah°f>k  +  Cp,k(T-T0) 

(13) 

Here,  A  is  the  area,  p  is  the  pressure,  H 
is  the  total  enthalpy  per  unit  mass,  Y*  is  the  mass 
fraction  of  the  kth  species  and  a;*  is  the  production 
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rate  of  the  kth  species.  The  right  hand  side  terms 
in  equation  5  represent  the  source  terms  due  to  area 
variation,  G^,  and  mass  addition  Gm-  me  is  the 
mass  flow  rate  of  the  injected  mass,  ue  is  the  velocity 
of  the  injected  mass,  He  is  the  total  enthalpy  per 
unit  mass  of  the  injected  mass  and  Yk>e  is  the  mass 
fractions  of  the  kth  species. 

The  governing  equations  described  here  are 
solved  by  a  simple  finite  difference  technique  using 
the  explicit  second  order  MacCormack  Scheme.  The 
solution  is  marched  time  accurately  to  steady  state. 

In  deriving  the  equations  2  and  3,  it  is  as¬ 
sumed  that  all  the  scales  of  motion  ranging  from  the 
smallest  Kolmogorov  scales  to  the  largest  integral 
length  scales  occur  in  the  flow[12]-[18].  This  implies 
that  in  a  numerical  simulation,  the  entire  range  of 
scales  should  be  completely  resolved.  However,  it 
has  been  found  that  this  is  not  always  necessary. 
Resolving  scales  larger  than  approximately  IO77  has 
been  found  to  yield  sufficiently  accurate  results[12], 
[15].  Hence,  in  the  present  work  also  only  scales 
larger  than  10 rj  are  resolved. 

4  Results 

A  schematic  of  the  dump  combustor  along  with  the 
close  up  of  the  injection  region  is  shown  in  figure  1. 
Figure  2  shows  the  schematic  of  the  model  that  rep¬ 
resents  the  dump.  Only  the  region  downstream  of 
the  dump  plane  is  modeled.  The  sudden  expansion 
is  modeled  as  a  section  with  a  large  increase  in  area. 
The  dump  area  then  changes  again,  the  diameter 
dropping  to  a  smaller  pipe.  This  causes  the  forma¬ 
tion  of  a  large  recirculation  area  behind  the  dump 
plane.  This  reduction  in  area  to  a  smaller  pipe  is 
modeled  as  a  straightening  of  the  one  dimensional 
channel  with  a  corresponding  decrease  in  cross  sec¬ 
tional  area.  This  kind  of  modeling  of  the  flow  is 
appropriate  because  we  are  only  trying  to  approxi¬ 
mate  the  mean  flow  field.  The  actual  dump  combus¬ 
tor  flow  field  contains  an  unsteady  shear  layer  which 
is  not  of  interest  here. 

In  order  to  determine  a  starting  point  for 
the  simulations,  the  data  from  an  LES  of  the  com¬ 
bustor  flow  field  and  the  corresponding  experimen¬ 
tal  investigation  were  used.  Figures  3  and  4  show 
the  experimentally  measured  mean  velociy  and  rms 
velocity  fluctuation  profiles  in  the  near  field  of  the 
dump.  The  center  line  of  the  shear  layer  is  identi¬ 
fied  as  the  line  of  maximum  rms  velocity.  A  plot 
of  the  mean  velocity  along  this  center  line  is  shown 
in  figure  5.  Also  shown  in  this  figure  is  the  velocity 
profile  obtained  by  varying  the  area  function  in  the 


model.  This  mean  axial  velocity  profile  is  used  in  all 
the  simulations. 

Next,  the  inputs  needed  for  the  stirring 
model  need  to  be  determined.  From  equations  2  and 
3,  it  is  clear  that  the  integral  length  scale(L)  varia¬ 
tion  needs  to  be  known  apriori.  In  ord<^Jto  calculate 
this  we  need  to  adopt  a  visualisation  of  the  shear 
layer  behind  the  dump.  The  picture  of  the  shear 
layer  adopted  in  this  work  is  shown  in  figure  6.  In 
the  near  field  of  the  dump,  the  shear  layer  is  assumed 
to  grow  linearly,  similar  to  a  mixing  layer  behind  a 
splitter  plate.  The  initial  layer  thickness  and  the 
linear  growth  rate  of  the  shear  layer  are  extracted 
from  the  LES  calculations.  The  LES  calculations 
show  that  in  the  near  field,  the  shear  layer  grows 
with  a  slope  of  about  0.2.  This  is  also  confirmed  by 
the  experimental  velocity  profiles  3,4.  This  thickness 
and  growth  rate  is  based  on  the  vorticity  thickness. 
Since  the  stirring  events  represent  the  action  of  ed¬ 
dies  of  different  sizes,  and,  the  vorticity  thickness 
represents  the  maximum  size  of  the  vortices  found 
in  the  shear  layer,  this  value  is  used  for  the  integral 
length  scale  L . 

Further  downstream  in  a  dump  combustor, 
the  shear  layer  curves  towards  the  wall  and  reat¬ 
taches  after  a  certain  length.  In  this  region  the  linear 
growth  rate  hypothesis  is  no  longer  accurate.  How¬ 
ever,  our  experience  with  this  geometry  and  the  ex¬ 
perimental  results  have  shown  that  the  primary  ac¬ 
tion  zone  is  in  the  near  field  of  the  dump  plane  and, 
hence,  it  is  expected  that  this  approximation  will 
not  significantly  alter  the  results.  Reasonably  good 
agreement  with  experimental  measurements  seem  to 
show  that  this  is  the  case. 

Another  parameter  that  needs  to  be  spec¬ 
ified  is  the  turbulent  Reynolds  number.  This  can 
be  done  by  specifying  the  turbulent  velocity  scale 
variation  in  the  shear  layer.  Since  the  shear  layer 
center  line  has  been  assumed  to  lie  along  the  line 
of  maximum  {7rm, ,  it  would  be  appropriate  to  spec¬ 
ify  the  Urms  as  the  turbulent  velocity  scale  in  the 
shear  layer.  From  figure  4,  we  see  that  this  maxi¬ 
mum  value  is  nearly  a  constant  along  the  length  of 
the  shear  layer  and  is  around  20%of  the  mean  ve¬ 
locity.  Hence,  in  the  present  work,  the  turbulent  ve¬ 
locity  is  specified  as  a  constant  20%  of  the  reference 
velocity. 

The  above  model  formulation  was  applied 
to  the  case  of  combustion  of  waste  surrogate  (Ben¬ 
zene)  in  a  dump  combustor.  One  of  the  primary 
quantities  of  interest  in  this  case  is  the  consumption 
efficiency,  also  known  as  the  Destruction  and  Re¬ 
moval  Efficiency  (DRE).  This  is  a  measure  of  how 
quickly  the  injected  waste  surrogate  is  consumed  in 
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the  dump.  A  quantitative  measure  of  this  is  given 
by  calculating  number  of  nines  as  follows, 


jjzof$  s  —  —logio(l.O  — 


%consumption 

loo 


(14) 


Hence,  a  consumption  of  99.99%  yields  a  DRE  of  4. 

A  plot  of  the  DRE  as  a  function  of  the  non- 
dimensional  axial  distance  is  shown  in  figure  7.  It 
is  seen  that  the  predicted  DREs  are  in  reasonable 
agreement  with  the  experimental  data  in  the  near 
field.  Further  downstream  the  DRE’s  are  overpre¬ 
dicted  by  the  model.  This  is  possibly  due  to  vari¬ 
ous  reasons.  The  chemical  mechanism  was  reduced 
from  a  larger,  more  detailed  mechanism  and  allows 
for  only  the  breakdown  of  the  benzene  molecules. 
There  is  now  evidence  to  suggest  that  benzene  may 
be  formed  as  an  intermediary  in  the  combustion  of 
some  hydrocarbons.  This  feature  is  not  included 
here. 

Another  reason  for  this  overprediction  is  the 
assumption  of  calorically  perfect  gas  (used  here  only 
as  an  approximation,  but  can  easily  be  relaxed)  used 
in  the  caloric  equation  of  state,  i.e.,  the  specific  heats 
of  the  gases  are  assumed  to  be  constant.  These  val¬ 
ues  are  assumed  to  remain  constant  at  their  values 
for  a  real  gas  at  300K.  However,  it  is  known  that  the 
specific  heats  can  change  by  a  factor  of  almost  3-3.5 
in  the  temperature  ranges  witnessed  in  the  present 
case  8.  This  results  in  a  overprediction  of  the  tem¬ 
perature  distribution  and  hence,  results  in  increased 
reaction  rates  and  higher  consumption  rates. 

Another  reason  for  this  increased  reaction 
rates  downstream  may  be  due  to  the  manner  in 
which  the  rearrangement  events  that  mimic  turbu¬ 
lent  stirring  are  implemented.  When  a  certain  loca¬ 
tion  is  chosen  for  the  stirring  event  and  the  size  of 
the  stirring  event  is  determined,  it  is  assumed  that 
equal  masses  are  stirred  and  that  these  masses  carry 
their  temperatures  with  them.  This  is  done  to  mimic 
the  increased  thermal  diffusivity  in  turbulent  flows. 
However,  in  the  one-dimensional  simulations,  since 
there  is  no  tranverse  thermal  diffusion  this  might 
result  in  a  slight  overprediction  of  the  thermal  dif¬ 
fusion  in  the  axial  direction.  This  results  in  higher 
temperatures  in  the  downstream  regions  and,  hence, 
higher  DRE’s. 

A  very  important  aspect  of  hydrocarbon 
combustion  is  the  emission  of  pollutants.  Efforts 
are  underway  to  reduce  pollution  by  the  emission 
of  dangerous  combustion  by-products  like  carbon- 
monoxide(CO)  and  oxides  of  nitrogen(ATOa;).  In  our 
Study,  one  of  the  major  criteria  used  to  validate  the 
model  has  been  the  prediction  of  pollutant  emis¬ 
sion.  A  plot  of  the  carbonmonoxide  mass  fraction 


(in  ppm)  is  shown  in  figure  9.  It  is  seen  that  for 
an  equivalence  ratio  of  0.54  (for  which  experimen¬ 
tal  data  was  available),  the  model  predictions  com¬ 
pare  very  well  with  the  experimental  data.  Also,  not 
much  variation  is  seen  in  the  near  field  with  change 
in  the  equivalence  ratio.  However,  further^down- 
stream,  significant  variation  is  seen. 

A  plot  of  the  variation  of  the  CO  emission 
with  change  in  the  equivalence  ratios  is  shown  in 
figure  10.  Values  plotted  are  the  CO  mass  frac¬ 
tion  (again  in  ppm)  at  X/D  =  3.0.  The  model 
predictions  are  in  excellent  agreement  with  the  ex¬ 
perimental  results.  The  CO  emission  is  seen  to  be 
minimal  around  equivalence  ratios  around  0.85:  At 
lower  and  higher  ratios,  the  emission  is  seen  to  rise 
significantly.  This  has  experimentally  been  observed 
to  be  caused  by  a  marked  decrease  in  flame  stability 
and  increase  in  its  sootiness  [19] 

Figure  11  shows  the  variation  of  NO  along 
the  non-dimensionalised  axial  coordinate.  The  NO 
production  mechanism  used  in  these  studies  is  the 
reduced  Zeldovich  mechanism  given  by  Chen  and 
Kollmann[9].  It  is  seen  that  the  agreement  between 
the  model  predictions  and  the  experimental  data  is 
reasonably  good.  However,  the  model  slightly  over¬ 
predicts  the  NO  emission.  This  is  due  to  the  fact 
that  the*  temperatures  are  slightly  over-predicted  as 
discussed  above.  This  further  indicates  that  the 
thermally  perfect  gas  model  for  the  caloric  equation 
of  state  needs  to  be  used.  The  CO  and  NO  emission 
are  seen  to  follow  very  different  trends.  While  CO 
decreases  with  downstream  distance,  NO  is  seen  to 
increase  with  X/D.  Figure  12  shows  the  variation 
in  NO  emission  with  equivalence  ratio.  Again,  in 
contrast  to  the  CO  variation  discussed  earlier,  NO 
emission  is  seen  to  increase  with  increase  in  <j).  These 
trends  are  captured  very  accurately  by  the  model. 

These  results  presented  here  show  that  the 
baseline  model,  as  presented  here,  is  capable  of  cap¬ 
turing  the  major  DRE  and  emission  trends  observed 
in  the  experimental  investigation.  However,  the 
global  research  effort  is  directed  towards  the  devel¬ 
opment  of  an  efficient  combustor.  Results  from  ex¬ 
periments  and  LES  investigations  show  that  acous¬ 
tic  and  mechanical  forcing  of  the  combustor  flow 
can  significantly  increase  the  combustion  efficiency 
in  the  combustor[10],[19].  This  forcing  has  also  been 
observed  to  reduce  the  CO  and  NO  emission  signif¬ 
icantly.  Our  future  work  will  be  directed  towards 
incorporating  a  mechanism  to  capture  these  effects 
into  the  model  presented  here. 
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5  Conclusion 

The  engineering  model  developed  is  seen  to  be  capa¬ 
ble  of  capturing  all  the  correct  trends  with  respect 
to  waste  consumption  and  pollutant  formations-  In 
most  cases,  not  only  are  the  qualitative  agreements 
good,  but  the  results  compare  quantitatively  well 
with  experimental  data.  The  model  is  fairly  simple 
to  apply  to  almost  any  kind  of  combustor  geometry 
and  has  not  required  any  significant  calibration  or 
tuning-  The  inputs  to  the  model  are  logically  ex¬ 
tracted  form  simple  experimental  and  LES  results. 
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Fig.  1.  Schematic  of  the  dump  combustor  used 
in  the  experiments  [19]  and  LES. 
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Fig.  2.  Schematic  of  the  model  geometry  used 
in  the  present  work. 


Fig.  4.  Axial  rms  velocity  profiles  measured 
experimentally  in  the  dump  combustor [  19] 
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Fig.  5.  Mean  axial  velocity  along  the  shear 
layer  centerline 
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Fig.  3.  Mean  axial  velocity  profiles  measured 
experimentally  in  the  dump  combustor!  19]. 


Fig.  6.  Topology  of  the  shear  layer  used  in  the 
present  work. 
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variation  along  the  axial  co-ordinate  with  the 
experimental  data. 


Fig.  8.  Temperature  profiles  predicted  by  the 
model  using  the  calorically  perfect  gas  model. 


Fig.  9.  CO  variation  along  the  axial  coordinate 
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Fig.  10.  CO  mass  fraction  variation  with  (p. 
Measurements  at  X/D  =3.0  are  shown 
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Fig.  11.  NO  variation  along  the  axial 
coordinate. 
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Fig.  12.  NO  mass  fraction  variation  with  (p. 
Measurements  at  X/D  =3.0  are  shown. 
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